J3RU 

1  9  APR  2000 


Development  of  Numerical  Constitutive  Models  for  Carbon-Carbon  Composites 

F49620-94-l-0259,P00002 
John  Whitcomb 

Center  for  Mechanics  of  Composites 
Texas  A  &  M  University 
-Final  Report 

May  15,  1994  -  August  31,  1998 


Objective 

Predict  the  stiffiiess  and  strength  of  oxidation-resistant  carbon-carbon  subjected  to  thermal 
and  mechanical  loads  in  an  oxidizing  environment 

Status  of  effort 

A  three-dimensional  progressive  failure  analysis  of  carbon-carbon  composites  was 
developed.  To  facilitate  the  study,  various  mesh  generators  and  graphical  pre-  and  post¬ 
processing  tools  were  developed.  These  tools  were  used  to  study  carbon-carbon  during  cool¬ 
down  from  processing  temperature  and  subsequent  mechanical  loads.  This  grant  was  tightly 
integrated  with  the  AASERT  grant  F49620-93- 1-0471,  which  ended  last  year.  The  student 
supported  under  the  AASERT  grant,  Clint  Chapman,  has  graduated  (and  is  employed)  with  a 
Ph.D.  This  final  year  of  the  subject  grant  has  been  invested  in  documenting  our  efforts  in  the  open 
literature  and  assembling  documentation  for  the  tools  developed  by  Dr.  Chapman  so  they  could  be 
used  in  subsequent  projects. 

Accomplishments 

Two  of  the  biggest  challenges  in  the  analysis  of  textile  composites  is  developing  a  valid  3D 
finite  element  mesh  and  deriving  boundary  conditions  for  the  smallest  possible  analysis  region 
(since  the  computational  challenge  is  inherently  high).  We  can  now  obtain  meshes  and  boundary 
conditions  for  plain  and  8  harness  satin  weaves  with  a  wide  range  of  waviness,  tow  cross-section, 
and  mesh  refinement.  Only  a  few  parameters  have  to  be  specified... the  rest  is  automatic.  We  have 
also  enhanced  our  visualization  programs  and  associated  utilities  so  that  it  is  much  more 
convenient  to  examine  stress  and  strain  distributions,  deformed  geometries,  damage  distribution, 
and  differences  between  models. 

We  have  continued  to  document  our  efforts  in  the  literature.  The  results  of  our 
simulations  were  presented  at  ICCE/5  in  July,  1998.  A  journal  paper  describing  the  simulations  is 
almost  complete.  The  title  is  “Thermally  Induced  Damage  Initiation  and  Growth  in  Carbon- 
Carbon  Composites.”  Another  paper,  “Derivation  of  Boundary  Conditions  for  Micromechanics 
Analyses  of  Plain  and  Satin  Weave  Composites,”  has  already  been  submitted  to  the  Journal  of 
Composite  Materials.  The  new  techniques  developed  in  this  paper  greatly  simplify  the  challenge 
of  deriving  boundary  conditions.  Reports  from  previous  years  of  this  grant  include  coDies  of 


£TXC  QUALITY  IM2PECTBD  4 


20000710  119 


REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searcl 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  sugges 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  oi  Management  and  Budget,  Paperw 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPOR.  ~  ^  -  vs* 

Final  Technical  Report  14  May  94  -  31  Aug  98 

4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

(U)  DEVELOPMENT  OF  NUMERICAL  CONSTITUTIVE  MODELS  FOR  F49620-94- 1-0259 

CARBON-CARBON  COMPOSITES 

2302/BS 

6.  AUTHOR(S)  61102F 

JOHN  WHITCOMB 


I  3.  REPOR  . 


AFRL-SR-BL-TR-OO- 


4.  TITLE  AND  SUBTITLE 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

TEXAS  A&M  UNIVERSITY 

CENTER  FOR  MECHANICS  OF  COMPOSITES 

COLLEGE  STATION,  TX  77843 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
AEROSPACE  AND  MATERIALS  SCIENCES  DIRECTORATE 
801  N.  RANDOLPH  STREET,  ROOM  732 
ARLINGTON,  VA  22203-1977 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


I  11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  IS  UNLIMITED 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

A  three-dimensional  progressive  failure  analysis  of  carbon-carbon  composites  was  developed.  To  facilitate  the  study,  various 
mesh  generators  and  graphical  pre-  and  post-processing  tools  were  developed.  These  tools  were  used  to  study  carbon-carbon 
during  cool-down  from  processing  temperature  and  subsequent  mechanical  loads.  This  grant  was  tightly  integrated  with  the 
AASERT  grant  F49620-93- 1-0471,  which  ended  last  year.  The  student  supported  under  the  AASERT  grant,  Clint 
Chapman,  has  graduated  (and  is  employed)  with  a  Ph.D.  This  final  year  of  the  subject  grant  has  been  invested  in 
documenting  our  efforts  in  the  open  literature  and  assembling  documentation  for  the  tools  developed  by  Dr.  Champan  so 
they  could  be  used  in  subsequent  projects. 


|  14.  SUBJECT  TERMS 


115.  NUMBER  OF  PAGES 


116.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2*89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DI0R,  Oct  94 


papers  and  Clint’s  PhD  dissertation,  so  they  will  not  be  included  herein.  Attached  is  a  copy  of  the 
paper  “Derivation  of  Boundary  Conditions  for  Micromechanics  Analyses  of  Plain  and  Satin 
Weave  Composites,”  which  has  been  submitted.  When  other  papers  are  ready  for  submission  to  a 
journal,  they  will  be  forwarded  to  the  technical  monitor. 


Peer  Reviewed  Publications 

1.  Whitcomb,  J.D.;  Chapman,  C.:  Effect  of  Assumed  Tow  Architecture  and  Mesh  Refinement 
on  the  Predicted  Moduli  and  Stresses  in  Plain  Weave  Composites.  Journal  of  Composite 
Materials.  Vol.  29,  No.  16,  pp.  2134-2159,  1995. 

2.  Whitcomb,  J.D.;  Srirengan,  K.;  Chapman,  C.:  Evaluation  of  Homogenization  for 
Global/Local  Stress  Analysis  of  Textile  Composites.  Composite  Structures.  Vol.  31,  No.  2, 
pp.  137-149,  1995. 

3.  Whitcomb,  J.D.;  Chapman,  C.:  Analysis  of  Plain  Weave  Composites  Subjected  to  Flexure. 
Mechanics  of  Composite  Materials  and  Structures.  Vol.  5,  No.  1,  pp.  41-53,  1998. 

4.  Srirengan,  K.;  Whitcomb,  J.D.;  Chapman,  C.:  Modal  Technique  for  Three-Dimensional 
Stress  Analysis  of  Plain  Weave  Composites.  Composite  Structures.  Vol.  39,  NO.  1-2,  pp.  145- 
156,  1997. 

5.  Whitcomb,  J.D.;  Chapman,  C.D.,  Tang,  X.:  Derivation  of  Boundary  Conditions  for 
Micromechanics  Analyses  of  Plain  and  Satin  Weave  Composites.  Submitted  8/98  to  Journal 
of  Composite  Materials. 

6.  Chapman,  C.D.;  Whitcomb,  J.D.:  Thermally  Induced  Damage  Initiation  and  Growth  in 
Carbon-Carbon  Composites  (manuscript  in  preparation). 


Conference  and  Workshop  Presentations 

1.  Whitcomb,  J.D.;  Srirengan,  K.;  and  Chapman,  C.:  Evaluation  of  Homogenization  for 
Global/Local  Stress  Analysis  of  Textile  Composites.  Presented  at  the 
AIAA/ASME/ASCE/AHS/ASC  35th  Structures,  Structural  Dynamics,  and  Materials 
conference,  Hilton  Head,  South  Carolina,  April  18-20,  1994. 

2.  Whitcomb,  J.D.;  Srirengan,  K.;  and  Chapman,  C.:  Simulation  of  Progressive  Failure  in  Plain 
Weave  Textile  Composites.  Presented  at  1994  ASME  Winter  Annual  Meeting,  Chicago, 
Illinois,  November  13-18,  1994. 

3.  Chapman,  C.;  Whitcomb,  J.D.;  and  Srirengan,  K.:  Analysis  of  Woven  Composites  Subjected 
to  Flexure.  Society  of  Engineering  Science  31*  Annual  Technical  Meeting,  College  Station, 
Texas,  October  10-12,  1994. 


4.  Srirengan,  K.;  Whitcomb,  J.D.;  and  Chapman,  C.:  Three-Dimensional  Analysis  of  Woven 
Composites  Subjected  to  Flexure.  Third  U.S.  National  Congress  of  Computational 
Mechanics,  June  12-14,  1995. 

5.  Chapman,  C.;  and  Whitcomb,  J.D.:  Thermomechanical  Characterization  of  Plain  and  Satin 
Weave  Carbon-Carbon  Composites.  Third  U.S.  National  Congress  of  Computational 
Mechanics,  June  12-14,  1995. 

6.  Whitcomb,  J.D.;  Srirengan,  K;  and  Chapman,  C.:  Modal  Technique  for  Three  Dimensional 
Stress  Analysis  of  Plain  Weave  Composites.  ICCM-10,  Whistler,  Canada,  August  14-29. 

1995. 

7.  Chapman,  C.;  and  Whitcomb,  J.D.:  Strategy  for  Modeling  Eight-Harness  Satin  Weave 
Carbon/Carbon  Composites  Subjected  to  Thermal  Loads.  Presented  at  37th 
AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural  Dynamics,  and  Materials  Conference, 
Salt  Lake  City,  UT,  April  15-17,  1996. 

8.  Whitcomb,  J.D.;  and  Chapman,  C..:  Thermally  Induced  Damage  Initiation  and  Growth  in 
Carbon-Carbon  Composites.  Presented  at  the  38th  AIAA/ASME/ASCE/AHS/ASC  Structures, 
Structural  Dynamics,  and  Materials  Conference,  Kissimmee,  FL,  April  7-10,  1997. 

9.  Whitcomb,  J.D.;  and  Chapman,  C.:  Damage  Initiation  and  Growth  in  Carbon-Carbon 
Composites  Subjected  to  Thermal  and  Mechanical  Loads.  Presented  at  the  1997  Joint 
ASME/ASEE/SES  Summer  Meeting,  Chicago,  IL,  June  28-  July  2,  1997. 

10.  Whitcomb,  J.D.;  Tang,  X.;  and  Chapman,  C.:  Simulation  of  Progressive  Failure  of  Woven 
Composites  Engineering,  to  be  presented  at  ICCE/5  Fifth  International  Conference  on 
Composites  Engineering,  Las  Vegas,  NV,  July  5-11,  1998. 

11.  Whitcomb,  J.D.;  Tang,  X.;  and  Chapman,  C.:  Comparison  of  Stress  Concentrations  in  Plain 
and  Satin  Weave  Composites,  to  be  presented  at  the  International  Mechanical  Engineering 
Congress  and  Exposition,  Dallas,  TX,  Nov.  1998. 


8/19/98  -  C:\whit_Rest\Papers\BC_forWeave\weaveBC_l.doc 


1 


Derivation  of  Boundary  Conditions  for  Micromechanics  Analyses  of  Plain  and  Satin  Weave 

Composites 
John  D.  Whitcomb 
Clinton  D.  Chapman 
Xiaodong  Tang 


Abstract 

Efficient  3D  analysis  of  periodic  structures  depends  on  identifying  the  smallest  region  to 
be  modeled  and  the  appropriate  boundary  conditions.  This  paper  describes  systematic  procedures 
for  deriving  the  boundary  conditions  for  general  periodic  structures.  These  procedures  are  then 
used  to  derive  the  boundary  conditions  for  plain  and  satin  weave  composites. 

Introduction 

Composite  materials  consist  of  a  combination  of  materials,  often  fibers  and  matrix.  The 
basic  challenge  of  micromechanics  is  to  determine  how  the  properties  and  spatial  distribution  of 
the  constituents  affect  the  overall  material  response.  The  overall  response  is  often  characterized 
in  terms  of  the  effective  engineering  properties,  such  as  extensional  moduli  and  Poisson's  ratio. 
These  are  also  referred  to  as  the  homogenized  engineering  properties.  To  obtain  these  properties 
a  representative  volume  element  (RVE)  or  unit  cell  is  identified  that  includes  all  characteristics 
of  the  composite.  The  literature  includes  a  wide  variety  of  analyses  for  various  unit  cells.  For 
example,  some  of  the  recent  numerical  studies  have  focused  on  hexagonal  and  square  arrays  of 
fibers  in  matrix  [1-4],  spherical  inclusions  in  matrix  [5],  and  textile  composites  [6-13],  The  latter 
configuration,  textile  composites,  presents  a  severe  challenge  to  the  analyst  because  of  the 
geometric  complexity.  Although  a  variety  of  analyses  exist,  including  some  very  detailed  three- 
dimensional  finite  element  analyses  [7-13],  there  has  not  been  a  systematic  description  of  how 
one  obtains  the  boundary  conditions  to  perform  the  analysis  of  textile  composites.  Reference  [14] 
gives  an  excellent  discussion  of  exploiting  symmetry  in  general.  However,  the  reference  does 
not  consider  periodicity  explicitly  or  composites.  Accordingly,  the  goal  of  this  paper  is  to  present 
a  systematic  procedure  for  deriving  the  boundary  conditions  for  both  full  unit  cell  and  partial  unit 
cell  analyses  of  plain  and  satin  weave  composites. 
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Background 

Implicit  in  micromechanics  analyses  is  that  a  small  region  of  the  microstructure  folly 
represents  the  behavior  of  a  much  larger  region  (usually  an  infinite  domain).  There  is  no  standard 
definition  of  this  region  in  the  literature.  Herein,  the  term  "unit  cell"  is  defined  to  mean  the 
smallest  region  that  represents  the  behavior  of  the  larger  region  without  any  mirroring  or  rotation 
transformations.  For  such  a  region  to  exist  there  must  be  a  basic  pattern  (of  geometry  and 
response,  such  as  strains,  stresses,  etc.)  that  is  repeated  periodically  throughout  a  domain.  In 
particular,  the  larger  region  can  be  synthesized  by  replication  and  translation  of  the  unit  cell  in 
the  three  coordinate  directions.  Figure  1  shows  examples  of  periodic  microstructure.  In  each 
case,  the  region  is  built  from  a  single  building  block  or  unit  cell.  Note  that  there  can  be  more  than 
one  definition  of  the  unit  cell  for  a  single  microstructure  (e.g.  case  (a)  in  Figure  1). 

To  determine  effective  engineering  properties,  the  periodic  array  is  subjected  to  a  series 
of  loads  consisting  of  either  macroscopically  constant  strain  or  stress.  The  term  macroscopically 
constant  indicates  that  the  volume  averaged  strain  for  every  unit  cell  is  identical.  (The  same  can 
be  said  for  the  stresses.)  Although  one  cannot  analyze  an  infinite  number  of  unit  cells,  it  is 
possible  to  derive  boundary  conditions  for  a  single  unit  cell  that  will  make  it  behave  as  though  it 
is  buried  within  an  infinite  array.  How  this  is  done  will  be  described  herein. 

It  is  convenient  to  think  of  the  challenge  as  consisting  of  two  parts.  The  first  part  consists 
of  identifying  the  appropriate  boundary  conditions  for  a  single  foil  unit  cell.  The  second 
challenge  is  to  exploit  the  symmetries  in  the  microstructure  of  the  unit  cell  and  the  loading  to 
reduce  the  analysis  region  to  just  a  fraction  of  the  foil  cell  in  order  to  reduce  the  computational 
costs. 

Derivation  of  basic  equations 

The  derivation  of  the  boundary  conditions  is  based  on  two  conditions.  The  first  is 
periodicity  and  the  second  is  equivalence  of  coordinate  system.  It  should  be  noted  that  that  these 
concepts  can  be  combined,  but  the  authors  choose  not  to  do  so.  The  periodic  conditions  state 
that  the  displacements  in  the  various  unit  cells  differ  only  by  constant  offsets  that-  depend  on  the 
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volume  averaged  displacement  gradients.  Further,  the  strains  and  stresses  are  identical  in  all  of 
the  unit  cells.  This  can  be  expressed  as 

«,(*«  +da)  =  Ui{Xa)  +  (j“^ 

Sii(xa+da)  =  eiJ{xa)  (1) 

o-v(xa+da)  =  crv{xa) 

where  dp  is  a  vector  of  periodicity.  This  vector  is  a  vector  from  a  point  in  one  unit  cell  to  an 
equivalent  point  in  another  unit  cell. 

These  conditions  are  sufficient  for  deriving  the  boundary  conditions  for  the  full  unit  cell. 
However,  it  is  usually  desirable  to  exploit  symmetries  in  order  to  be  able  to  analyze  a  smaller 
region.  The  concept  of  "Equivalent  Coordinate  Systems"  (ECS)  is  useful  in  identifying  the 
symmetries  and  constraint  conditions.  Coordinate  systems  are  equivalent  if  the  geometry,  spatial 
distribution  of  material,  loading,  and  the  various  fields  that  describe  the  response  (e.g. 
displacement,  strains,  etc.)  are  identical  in  the  two  systems. 

For  example,  x{  and  xt  are  equivalent  coordinate  systems  if 

«/(*«)  =  «,  CO 
*#(*«)  =  **(*«) 

*>(*«)  =  ^(*«)  (2) 


etc. 

A  visual  technique  for  determining  whether  a  new  coordinate  system  is  equivalent  is  to  draw  it 
on  a  three-dimensional  view  of  the  body.  This  view  should  include  the  load  vectors  and  a 
description  of  the  spatial  variation  of  the  material  properties.  If  the  body  can  be  rotated  (i.e.  the 
view  changed)  such  that  the  new  view  and  new  coordinate  system  look  identical  to  the  original 
view  and  coordinate  system,  the  coordinate  system  is  equivalent. 

Exploiting  equivalent  coordinate  systems  requires  identification  of  coordinate  transformations 
that  provide  useful  constraint  conditions.  Consider  the  equivalent  coordinate  system  xjas  a 
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potentially  useful  ECS.  Define:  =  direction  cosines  for  transformation  from  the  original  (xy) 
to  the  equivalent  (xj)  coordinate  system.... i.e.  x,  =  atj  xy 

=  a(iUj{Xo) 

=  (3) 

^lJ(xa)  =  aimajnam„(xa) 

Combine  equations  (2)  and  (3)  to  give 

",  CO  =  <¥',(*«) 

=  (4) 

Finally,  replace  xa  with  uray  xt  so  that  everything  is  now  expressed  in  terms  of  a  single 
coordinate  system. 

»,(*«)  s<V'/(a**y) 

(5) 


Sometimes  it  is  necessary  to  switch  the  sign  of  all  of  the  loads  to  make  the  two  coordinate 
systems  equivalent.  If  tension  and  compression  properties  are  different,  one  must  not  switch  the 
sign  of  the  load  to  obtain  an  equivalent  coordinate  system. 


To  generalize  equations  (5),  a  factor  y  is  introduced. 

y- 1  if  load  reversal  is  not  required 
y  =  - 1  if  load  reversal  is  required 
Equation  (5)  then  becomes 

«i(xa)  =  yaiiuJ(a<9.xJ) 

ey(Xa)  “  Yaimaine«m{a<f  Xj) 
^M  =  yaimajn(Tm„(acg.xJ) 


K 


K 


(6) 
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The  constraints  that  are  implied  by  equations  (6)  are  very  useful  in  deriving  the  required 
boundary  conditions.  This  will  be  demonstrated  later  through  examples. 

The  value  of  the  factor  y  is  easily  determined  from  the  loading,  since  only  simple  loading 
is  considered.  In  particular,  the  model  is  subjected  to  either  macroscopically  constant  strain  or 
stress.  For  example,  one  might  specify  a  volume-averaged  value  of<r12  (with  the  other  volume- 
averaged  stresses=0).  The  sign  of  y  can  be  obtained  by  examining  the  effect  of  the  coordinate 
transformation  on  the  sign  of  crl2 .  In  fact,  we  can  consider  all  six  load  cases  at  once.  For 
example,  if  the  coordinate  transformation  is  mirroring  about  the  x,  =0  plane. 


■-1 

0 

0* 

"  CTii 

~CT12 

1 - 

fO 

b 

I 

aij  - 

0 

1 

0 

which  gives  = 

-<*12 

CT22 

°23 

0 

0 

1_ 

ra .3 

^23 

*33  . 

Accordingly,  the  values  of  y  are  as  follows  for  the  six  different  load  conditions  for  mirroring 
about  the  x,  =0  plane. 

Stress  Cjj  cr22  &33  o^2  Cjj  cr13 
y  111-11-1 

The  values  of  y  for  six  different  coordinate  transformations  and  the  six  load  cases  are 
summarized  in  Table  1  below.  The  values  of  y  for  specified  strains  e9  or  displacement  gradients 

-^-are  the  same  as  for  stress  loading. 

dXj 


Table  1  Values  of  y  for  various  transformations  and  load  cases. 


Transformation  Load  Case 


cr22 

^33 

^12 

^23 

^13 

Mirror  about  xl  =  0 

1 

1 

1 

-1 

1 

-1 

Mirror  about  x2  =  0 

1 

1 

1 

-1 

-1 

1 

Mirror  about  x3  =  0 

1 

1 

1 

1 

-1 

-1 

^•rotation  about  xl 

1 

1 

1 

-1 

1 

-1 

^rotation  about  x2 

1 

1 

1 

-1 

-1 

1 

/r  rotation  about  x3 

1 

1 

1 

1 

-1 

-1 

6 


8/19/98  -  C:\whit_Rest\Papers\BC_forWeave\weaveBC_l.doc 


The  values  of  y  also  indicate  what  types  of  loading  can  exist  for  a  particular  type  of  equivalent 
coordinate  system.  For  example,  for  mirroring  about  the  x,  axis  one  cannot  impose  non-zero 
volume  averaged  values  of  both  sxx  and  en  at  the  same  time  ...  the  symmetries  are 
incompatible  since  they  require  different  values  of  y .  Furthermore,  if  one  imposes  a  non¬ 
zero  (f,,)  and  retains  the  other  strains  as  unknowns  to  be  determined  as  part  of  the  solution, 
and  (£,3)  must  be  found  to  be  zero. 

Boundary  conditions  for  symmetrically  stacked  plain  weave 

If  one  is  willing  to  analyze  the  entire  unit  cell,  the  periodic  conditions  for  the 
displacements  can  be  used  to  specify  the  boundary  conditions.  However,  to  reduce  the  region  to 
be  analyzed  requires  a  bit  more  effort.  The  following  describes  a  systematic  way  to  go  about 
accomplishing  this  task.  Although  the  focus  is  on  the  plain  weave,  the  techniques  are  quite 
general.  Figure  2a  is  a  unit  cell  for  a  symmetrically  stacked  plain  weave  composite.  As  can  be 
seen  from  the  figure,  the  planes  xt  =  0  are  planes  of  geometric  symmetry. 


The  first  step  will  be  to  develop  relationships  for  the  planes  xx  =  ±a .  (vector  of  periodicity 
da  =  [2a,0,0])  Substituting  xx  =  -a  into  the  periodic  relations  gives 


Ui(a,  x2,  x3)  =  x2,x3)  +  \f~-/2<3 

\  COC^  j 


£ij(a>x2,x3)  =  £  if  (-a,  x2  ,x3) 
er,j(a,x2,x3)  =  aiJ(-a,x2,x3) 


(7) 


Consider  the  equivalent  coordinate  system  (ECS)  xa  that  is  obtained  by  mirroring  about  xx  =  0 . 
The  direction  cosines  for  this  transformation  are 


av  = 


-1  0  O' 
0  1  0 
0  0  1 


(8) 


Substituting  the  ai}  above  into  equations  (6)  yields  the  following  relationships  for  displacements 
and  stresses. 
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"i  -M, 

u2  =y  u2 

ui -Wjjr,)  L 

O^ll  °12  ®"l3  °ll  —0r12  ~^13 

o’  12  0"22  CTjj  =  —  <Tj2  0"22  ®23 

cr13  0-33  ^33  J(a^C2^3)  L_cr13  °23  ^33 


(9) 


Tractions  will  transform  like  the  displacements  and  strains  like  the  stresses,  so  the  details  are  not 
shown.  These  periodic  and  ECS  conditions  will  now  be  combined  ...  first  for  displacements  and 
then  for  the  stresses. 


Recall  that  the  value  of  y  depends  on  whether  the  loads  must  be  reversed  for  the  new  coordinate 
system  to  be  equivalent.  The  constraints  indicated  by  equation  (11)  are  listed  below.  The  ones 
that  provide  useful  information  are  boxed.  The  remaining  constraints  are  not  useful  for 
performing  the  analysis,  but  they  can  be  used  for  evaluating  the  results. 
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This  can  be  expressed  quite  succinctly  as  follows 

•  If  y  =  1 ,  normal  displacements  on  xx  =  -a  are  equal.  Because  of  periodicity,  the  normal 
displacements  on  x1  =  a  are  also  equal  and  opposite  those  onxj  =  -a . 

•  If  y  —  -1 ,  tangential  displacements  are  constant  in  each  direction  on  Xj  =  -a .  Because  of 
periodicity,  the  tangential  displacements  in  each  direction  on  xx  =  a  are  also  equal  and 
opposite  the  corresponding  displacements  onxl= -a. 

Stresses 

Now  consider  the  constraints  on  the  stresses.  Combine  equations  (7)  and  (9)  to  obtain  the 

following  useful  relations  (the  remainder  are  trivially  satisfied  or  do  not  affect  the  boundary 

conditions.) 

r  =  1:  o-12(±a,x2,x3)  =  0  crl3(±a,x2,x3)  =  0 

r  =  ~  1:  crn(±a,x2,x3)  =  0 

These  can  be  summarized  as  follows  (using  Cauchy's  stress  formula,  rt  =  afinj) 

•  If  /  =  1  shear  traction  on  xx  =  ±a  is  zero 
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•  If  y  =  —l  normal  traction  on  xt  =  ±a  is  zero 


Because  of  the  symmetries  in  the  microstructure,  we  will  obtain  analogous  results  for 
x2  =  ±a  and  x3  =±t .  That  is,  on  these  planes 

•  lfy  =  l,  normal  displacements  are  equal,  and  shear  tractions  are  zero. 

♦  Hy  =  -1,  tangential  displacements  are  equal,  and  normal  tractions  are  zero. 

Next  we  proceed  to  exploit  symmetries  within  the  cell  to  reduce  the  region  which  must  be 
modeled.  Again,  we  will  consider  an  ECS  obtained  by  mirroring  about  the  xl  =  0  plane. 
However,  this  time  we  examine  the  displacements,  etc.,  along  the  plane  =  0 .  Equations  (6) 
become 


O’n 

V J2 

^13  " 

’  0-lt 

-^12 

-^3 

^12 

^22 

C23 

=  r  -°X7 

°22 

Cr23 

3 

^23 

°"33. 

(O.xj.x,)  -  CTl3 

^23 

^33 

Therefore,  on  x3  =  0  we  have  the  following  constraints 

r= 1  r~- 1 

normal  displacement  =  0  tangential  displacements  =  0  (15) 

tangential  tractions  =  0  normal  tractions  =  0 

If  we  repeat  this  exercise  for  the  planes  x2=0  and  x3  =  0  (which  are  also  planes  of 
microstructural  symmetry),  we  will  obtain  analogous  results.  We  are  now  down  to  1/8  of  the 
original  unit  cell  (Figure  2(b)).  The  analysis  region  can  be  reduced  to  1/32  of  the  unit  cell.  To 
simplify  the  derivation  of  the  boundary  conditions  for  the  1/32  unit  cell,  the  coordinate  system 
origin  is  shifted  to  the  center  of  the  1/8  unit  cell  (Figure  2(c)). 

Consider  a  rotation  about  the  xl  axis  of  180°  to  obtain  an  ECS. 
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’10  0 
«„= 0  -i  0 
0  0-1 


The  periodic  conditions  in  equation  (6)  yield  (specialized  for  x2  —  0). 


= r\  ~u. 


3  JCx^Ojcj) 


°12  °"l3 

°"l2  °22  ^23 

°23  Or33j(,l.0.,J) 


3  J(*i,0.-*j) 


=  -cr, 


The  results  in  the  following  displacement  constraints 

_ Tzi _ 

^1(^1  A^)  A  "^3 ) 
w2  (*1 »°» *3)  =  “W2  (*1  A-*3  ) 

^3(^1  A  ^3)  -  “^(^iA1  ^3) 


^12  an 


°33  J^.O.-x,) 


_ y=-i _ 

«l(*lA*3)  =  -«l(*lA-*3) 
«2(XlAX3)  =  «2(X|A-^) 
W3(X1A^3)  =  W3(XlA-^3) 


Note  that  constraints  are  derived  for  all  three-displacement  components.  Hence,  traction 
boundary  conditions  do  not  have  to  be  considered.  Now  one  can  replace  the  region  x2  >  0  with 
these  constraints.  The  analysis  region  is  now  1/16  of  the  unit  cell. 

Now  consider  a  rotation  about  the  x2  axis  of  180°  to  obtain  an  ECS  and  specialize  equation  (6) 
for  the  plane  xx  =  0 . 


-1  0  0 

0  1  0 

0  0-1 


The  result  is 


=r\  u. 


3 -*(0^2  -*j) 
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O'u 

^13* 

'0-11 

-*!2 

0*13  ' 

°X2 

^22 

^23 

=r 

-O’u 

0*22 

“0-23 

.O’, 3 

0-23 

u» 

1 _ 

.0-13 

“0-23 

£ 

w 

(20) 


This  results  in  the  following  displacement  constraints 


_ r  =  i _ 

w,  (0,  x2 ,  x3 )  =  -w,  (0,  x2  -x3 ) 

^2  *2  >  *3 )  =  ^2  ^2 ^3  ) 

w3(0,X2,x3)  =  -w3(0,x2,-x3) 


y  =  - 1 _ 

Wl(0.^2»X3)=Wl(0>X2  -^3> 

u2(0,x2,x3)  =  -«2(0,x2 -x3) 
z/3  (0,  x2,x3)  =  w3  (0,  x2  ,-x3 ) 


(21) 


Again,  one  obtains  constraints  for  all  three-displacement  components.  Hence,  the  region  x,  >  0 
can  be  reduced  with  these  constraints,  leaving  an  analysis  region  of  only  1/32  of  the  unit  cell. 


Table  2  summarizes  the  boundary  conditions  for  symmetrically  and  simply  stacked  plain 
weaves  under  extension,  cr12  (or  en  )  shear,  and  cr13  (or  sl3 )  shear  load  conditions. 


Boundary  conditions  for  satin  weave 

The  boundary  conditions  for  a  simply  stacked  8-hamess  satin  weave  (Figure  3)  will  be 
derived  in  this  section.  Modifications  for  a  symmetrically  stacked  satin  weave  are  also  given.  If 
the  entire  unit  cell  is  to  be  analyzed,  the  periodic  conditions  given  in  equation  (1)  are  sufficient. 
The  vectors  of  periodicity  to  be  used  in  equation  (1)  are  [3a,-a,0],  [2a,2a,0],  [-a,3a,0],  [0,0,t]. 
However,  it  is  possible  to  analyze  half  of  the  unit  cell  if  both  the  periodicity  and  an  equivalent 
coordinate  system  are  exploited.  The  half-cell  that  will  eventually  be  used  is  shown  in  Figure 
3(b).  The  equivalent  coordinate  system  to  be  used  is  obtained  by  a  rotation  of  n  about  the  x3axis. 
For  such  a  transformation  equation  (6)  gives  the  following  constraints  for  the  full  unit  cell. 

«i(x,,x2,x3)  =  -^,(-xI,-x2,x3) 

m2(x,  ,  X2 ,  x3)  =  -yu2  (-x,  ,-x2 ,  x3)  (22) 

w3(x„x2,x3)  =  >«3(-x„-x2,x3) 

Boundary  conditions  for  the  plane  x2  =  0  can  be  obtained  immediately  by  substituting 
x2  =  0  into  equation  (22).  Equation  (22)  will  also  be  useful  later  as  each  pair  of  faces  of  the  unit 
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cell  is  considered.  Suitable  constraints  must  be  derived  for  the  other  faces  in  order  to  reduce  the 
analysis  region  to  one-half  of  the  unit  cell. 

Faces  x1  =  ±^- 
2 

The  vector  of  periodicity  that  connects  these  faces  is  da  =  [3a, -a,  0] .  Substitution  of  this 
vector  into  equation  (1)  gives 

"‘(t’*2 ~a’x*)= w{■T,X2,X3)+3a(^^}“a(^^},  7 = 1,2,3  (23 

cl  3a 
where  —  <x2<  — . 

2  2  2 

The  constraints  in  equation  (23)  are  used  to  specify  the  boundary  conditions  for  the  shaded 
regions  in  Figure  4(a). 


Combining  equation  (22)  and  (23),  gives  boundary  conditions  for  the  face  x,  =  — as 

2 


follows: 


(3a  ^  (3a  ^ ,  o  \  l^u , 

-°-x’  r FF^rW, 

(3a  \  (3a  \  /8m2\  /du 

^  J- J+3a\&7_a\£ 

(3 a  (3a  V,  / 0w3 \  /du3\ 


where  —  <x2< — .  Note  that  these  constraints  only  involve  points  on  the  plane  xl  =  — . 

2  2  2 

Equation  (24)  is  used  to  specify  the  boundary  conditions  for  the  points  specified  by  0^x2  <a  in 
equation  (24).  (shaded  region  in  Figure  4(b)  ). 


Faces  x,  +  x2  =  ±2a  (beveled  edges) 

The  vector  of  periodicity  that  connects  these  faces  is  da  =  [2a,2a,0] .  Substitution  of  this 
vector  into  equation  (1)  gives 


8/19/98  -  C:\wMt_Rest\Papcrs\BC_forWeave\weaveBC_l.doc 


13 


wi(  x2,x2+2afx3 )  —  **/(  x2  2a,  x2,x3) +2a^_^  ^+2a^_^  /  1,2,3 

(2f 

Combining  equation  (25)  and  the  equivalent  coordinate  system  (22)  gives  the  following 
boundary  conditions  for  face  x,  +x2  =  -2a . 

—yux{x2,—x2  —  2a,x3)  =  ul(—x2  —  2a,x2,x3)  +  2a(^~  -^+2a^—  ^ 

-^(*2 -*2 ~2a,x3)  =  u2(-x2-2a,x2,x3)+2a(^^+2a(^j  (26 

2^3 (^2 »  ■ *2  “ 2a, x3)  =  u3{—x2  —  2a,x2,x3)  +  2a^j— ^+2a^— ^ 

where  — —  <  x,  <  — .  The  constraints  in  equation  (26)  only  involve  points  on  the  plane 
2  2 

x,  +  x2  =  -2a .  This  can  be  verified  by  observing  that  the  coordinates  satisfy  the  condition 
x,  +  x2  =  -2a .  Similar  boundary  conditions  can  also  be  obtained  for  face  X;+X2=2a,  but  they  are 
not  needed  for  the  half  unit  cell. 


Faces  x2  =  ±: 


The  vector  of  periodicity  that  connects  these  faces  is  da  =[-a,3a,0] .  Substitution  of  this 
vector  into  equation  (1)  gives 

w{x,,T,X3)=M{Xl"T,X3)"a(^)+3a(^}’  /= 1,2,3  (27; 

Combining  equation  (27)  and  the  equivalent  coordinate  system  (22)  gives  the  following 

3a 

boundary  conditions  for  face  x2  = - 
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(  3a  ^  (  3a  }  /duA  /dux 

-W^x,+a,-T,x3J  =  ^x„-T,)r, 

(  3a  \  (  3a  /3m2  \  /ch/2 

^2[  X'+a’  2’X3JU\Xl’  2  ’*3J  AdxJ  3  \dx2 

(  3a  \  (  3a  V  /du3\  /^/3\ 

+°--T’x>n{x'-~-xTa^rat^) 


where  <*.<  —  .  Similar  boundary  conditions  can  alse  be  obtained  for  face  xf=3zJ2,  but 
2  2 

they  are  not  needed  herein. 


Faces  x,  =  ±— 
3  2 


The  vector  of  periodicity  that  connects  these  faces  is  da  -  [0,0,  /] .  Substitution  of  this 


vector  into  equation  (1)  gives 


<x„xii)  =  ,(x„x!,-i) 


+  t(^\  ^  =  1,2,3 


Up  to  this  point  the  boundary  conditions  described  have  been  valid  for  both  simple  and 
symmetric  stacking  except  equation  (29)  is  valid  for  simple  stacking  only.  For  symmetric 

stacking,  the  boundary  conditions  on  x3  =  ±—  are  the  same  as  for  the  symmetrically  stacked 

2 

plain  weave  discussed  earlier,  that  is 

•  If/  =  1 ,  normal  displacements  are  equal,  and  shear  tractions  are  zero. 

♦  If/  =  -1 ,  tangential  displacements  are  equal,  and  normal  tractions  are  zero. 

Comments  on  boundary  conditions 

Earlier  in  the  discussion  of  Table  1  it  was  mentioned  that  certain  macroscopic  modes  of 
deformation  are  incompatible  when  transformations  are  used  to  reduce  the  analysis  region  to  less 
than  the  full  unit  cell.  Although  all  the  macroscopic  displacement  gradients  were  listed  in  the 
boundary  conditions,  whenever  a  particular  load  case  is  considered,  many  of  the  gradients  must 
be  zero.  Also,  it  is  necessary  to  eliminate  rigid  body  motion  of  the  model.  Rigid  body 
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translations  are  prevented  by  constraining  a  single  point  in  all  three  coordinate  directions.  Rigid 
body  rotations  are  prevented  by  imposing 


for  i^j 


) 


This  causes  no  loss  of  generality  in  the  modes  of  deformation.  The  simplified  boundary 
conditions  for  simply  and  symmetrically  stacked  8-hamess  satin  weaves  are  summarized  in 
Table  3. 


As  listed,  specifying  the  magnitude  of  the  displacement  gradients  specifies  the  magnitude 
of  the  load.  If  one  specifies  one  displacement  gradient  to  be  non-zero  and  the  rest  to  be  zero,  then 
obviously  there  is  only  one  non-zero  volume  averaged  strain,  but  there  are  generally  several  non¬ 
zero  volume  averaged  stresses.  Conversely,  if  one  displacement  gradient  is  specified  and  the  rest 
are  determined  as  part  of  the  solution,  there  will  generally  be  several  volume  averaged  strains, 
but  only  one  volume  averaged  stress. 


Conclusions 

A  systematic  procedure  was  developed  for  deriving  boundary  conditions  for  partial  unit 
cell  models  of  materials  with  periodic  microstructure.  The  procedure  combines  the  concepts  of 
periodicity  and  equivalent  coordinate  systems.  Boundary  conditions  were  derived  for  simply  and 
symmetrically  stacked  plain  and  8-hamess  satin  weave  composites  for  various  types  of 
macroscopic  loads. 
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Table  2  Summary  of  boundary  conditions  for  plain 


Table  3  Summary  of  boundary  conditions  for  8-hamess  satin  weave 


Unit  cell 

(b)  Simply  stacked  satin  weave 


Figure  1.  Examples  of  periodic  microstructure 


(a)  Full  unit  cell 


(b)  1/8  unit  cell  considered  (shown 
cut  from  full  unit  cell) 


(c)  New  coordinate  system  (4)  1/32  unit  cell 


Figure  2  Coordinate  systems  and  regions  used  in  derivation  of 
boundary  conditions  for  symmetrically  stacked  plain  weave 


Figure  3  Coordinate  system  and  regions  used  in  derivation  of 
boundary  conditions  for  simply  stacked  8-harness  satin  weave 


(a)  Regions  where  boundary  conditions 
are  based  on  periodic  condtions. 


(b)  Region  where  boundary  conditions  are 
based  on  both  periodic  and  equivalent 
coordinate  system  constraints. 


Figure  4  Paired  regions  for  multipoint  constraints . 


